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Abstract 

Bogoliubov-de Gennes (BdG) equations and the excitation spectrum of a two- 
component Bose-Einstein condensate (BEC) are derived with an arbitrary in- 
teraction between bosons, including long-range and short range forces. The 
nonconverting BEC mixture segregates into two phases for some two-body 
interactions. Gross-Pitaevskii (GP) equations are solved for the phase segre- 
gated BEC. A possibility of boundary-surface and other localised excitations 
is studied. 
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Neutral and charged (Coulomb) Bose-gases became recently of particular interest motivated 
by the observations of BEC in an alkali vapor and by the bipolaron theory of high tem- 
perature superconductors 0, respectively. Their theoretical understanding is based on the 
Bogoliubov 1^ displacement transformation, separating a large matrix element of the conden- 
sate field operator, 0(r, t), from the total ■?/'(r, t) and treating the rest '?/'(r, t) = ■ip{r, t)— 0(r, t) 
as a small fluctuation. The resulting (Gross-Pitaevskii) equation for the condensate wave 
function 0(r, t) provides the mean-field description of the ground state and of the excitation 
spectrum |8HlO|. Beyond the mean- field approach the Bogoliubov-de Gennes (BdG)-type 
equations were derived [0, describing eigenstates of the 'supra'condensate bosons. 

In this letter we extend the Bogoliubov theory to the two-component nonconverting 
condensate by deriving GP and BdG equations and excitation spectrum with an arbitrary 
two-body interaction, and solving GP equations for a phase-segregated BEC. Our motivation 
originates in the recent experimental and theoretical P,p!0| studies of BEC of ^^Rb atoms 



in two different hyperfine states, and also in an observation ||T2[ that the condensation tem- 
perature of a two-component charged Bose gas quantitatively describes the superconducting 
critical temperature of many cuprates. Differently from Ref. we consider nonconvert- 

ing components (such as different elements), and an arbitrary (rather than short-ranged) 
two-body interaction. 

The Hamiltonian of the two-component (1 and 2) mixture of bosons in an external field 
with the vector, A{r,t), and scalar, Uj{r,t), potentials is given by 

(V-zg,A(r,t))2 



i=i,2 



2mj 

+ lY.JdrJ dv%,,{v - r')^](r)^,(r)4(r')^,v(r'), (1) 

where mj and qj are the mass, and the effective charge (if any) of the boson j {h = 1). If 
the two-body interactions ^/(r) are weak, the occupation numbers of one-particle states 
are not very much different from those in the ideal Bose-gas. In particular the lowest energy 
state remains to be macroscopically occupied and the corresponding component of the field 
operator ^/'(r) has anomalously large matrix element between the ground states of the system 



containing + 1 and bosons. Hence, it is convenient to consider a grand canonical 
ansamble, introducing the chemical potentials fij to deal with the anomalous averages (f){r) 
rather than with the ofF-diagonal matrix elements. Using the Bogoliubov displacement 
transformation in the equation of motion for the field operators '?/'(r, t) = 0(r, t) + ip{r,t) 
and collecting c — number terms of 0, and terms linear in the supracondensate boson 
operators ip one obtains a set of the GP and BdG-type equations 0,^|lOl. The macroscopic 
condensate wave functions (i.e the order parameters) obey two coupled GP equations 

2-0,(r,t) = /i,0,(r,t) + 5: 1 dv'V,Ar-r')\M^',t)\'^,ir,t) (2) 

f 

with the single-particle Hamiltonian hj = —(V — iqjA{r,t))'^ /2mj + Uj{r,t) — fij. The 
supracondensate wave-functions satisfy four BdG equations 

Y: I dr%,> (r - r') [| (r', t) \ \, (r, t) + (r', (r, t)u,> (r', t) + 

d 

(j)j{r',t)(f)j>{r,t)vj>{r',t)] = i—Uj{r,t) - hjUj{r,t), (3) 

and 

Y: J dv%,\v - r')[|0,v(r',t)|\-(r,t) + 0,(r', t)0*,(r, t)t;,,(r', t) + 

3' 

<P]{v\t)<P),{v,t)u,,{v\t)] = -i^^v,{v,t) - h*v,{v,t). (4) 
Here we have applied the linear Bogoliubov transformation for 

= I]«ni(r,t)(a„ + /3n) +<^-(r,t)(aJ, + /5^), (5) 

n 

where a„, j3n are bosonic operators annihilating quasiparticles in the quantum state n. There 
is a sum rule, 

^K(r,t)<^.(r',t) - t;„,(r,t)<^.(r',t)] = 5{v - r'), (6) 

n 

which retains the Bose commutation relations for all operators. 

In the homogeneous system with no external fields the excitation wave functions are 
plane waves 



Mkj(r, t) = exp[ik ■ r — iE(k)t], 



(7) 



t>kj(r, t) = exp[ik ■ r — iE{k)t], 



(8) 



The condensate wave function is {r,t) independent in this case, (j)j{r,t) = (pj. Solving two 
GP equations one obtains the chemical potentials as 



fii = Vrii + Wn2 
^2 = Un2 + Wrii, 

and solving four BdG equations one determines the excitation spectrum, -E(k) from 
V^<Pf 6(k) + ^(k) W^<Pl<Pl W^<Pl<P2 



(9) 
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Here ^i(k) = A;^/2mi + Vk'^i,^2(k) = k'^/2m2 + U\ji2, and Vk, f/k, Wk are the Fourier com- 
ponents of V^ii(r), V22(r) and ^42(1"), respectively, V = Vq, U = Uq, W = Wq, and Uj = \(f)j\'^ 

the condensate densities. There are two brunches of excitations with the dispersion 

1/2 



Ei,2(k) = el{k) + el{k) ± J[el{k) - ei(k)]2 + Wln,n2 



mim2 



(11) 



where ei(k) = \Jk^/{Ami) + fc^Vk^i/mi and e2(k) = ^k^/^Aml) + k'^U\ji2/m2 are Bogoli- 
ubov's modes of two components. If the interactions are short-ranged (hard-core) repulsions, 
so that Vk = V^, f/k = U, Wk = the spectrum, Eq.(ll) is that of Ref. In the long- 
wave limit both brunches are sound-like with the sound velocities 

/ I \ 1/2 

si,2 = 2"^/^ [yrii/mi + Un2/m2 ± \l\Vni/mi - Un2/m2Y + 4W^^nin2/(mim2) ) . (12) 



The lowest brunch becomes unstable (s2 < 0) if 1^ > yUV . When W = VUV this brunch 
is quadratic in the long-wave limit 

k^ ' 



^2(k) 



Vnirrii + Un2m2 
2{mim2Y/'^ V Vnim2 + Un2'mi 
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(13) 



If mi = 7712 = 771 it becomes 'collisionless' |]T0[, i-e -E'2(k) = A;^/2m. 

The finite-ranged interactions drastically change the whole spectrum. In the extreme 
case of the long-range Coulomb interaction, l^(k) = ATrql/k"^, ?7(k) = ATtq^/k'^., and iy(k) = 
4vrgig2/^^ the upper brunch is the geometric sum of the familiar plasmon modes for 
A; ^ 0, 



Ei(k) = Jl^ + l^, (14) 

V Till 77l2 



while the lowest brunch is 



qlnimi + ql7i277i2 

(15) 



2(mim2)^/^ ^ qinim2 + 

Remarkably, this mode is 'collisionless' at a7iy charges of the components if mi = m2, 
Fig.l. Physically, it describes neutral oscillations of condensates, when a charge fluctuation 
in one component is nullified by a charge fluctuation in another component, and the total 
charge density does not fluctuate. We conclude that while the Coulomb Bose condensate is 
a superfluid (according to the Landau criterion), mixture of two Coulomb Bose condensates 
is not. In a more general case the interaction might include both the long-range repulsion 
and the hard-core interaction as in the case of bipolarons or any other preformed bosonic 
pairs [^]. Combination of both interactions, i.e. Vk, W^k oc co7ista7it + l/fc^ transforms 
the lowest quadratic mode into the Bogoliubov sound. Hence, two component condensate 
of bipolarons is a superfluid. 

Finally, let us discuss the phase segregated state of the two-component nonconverting 
mixture with the hard-core interactions when W > VUV [|14|. In that case chemical poten- 
tials determined in Eq.(9) are no longer correct. Minimizing the free energy with respect to 



the equlibrium concentration we find the densities ti'^ = 7ii + ^JU/V7i2 and rig = 712 + yV/Uni 
of two separated phases. The phase boundary is described by two coupled one-dimesional 
CP equations, Eq.(2), 



r + h-fi-rfj^ = 0, (16) 
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and 



«^ + /2 - /I - rhf. = 0. (17) 



where we introduce two real dimensionless order parameters, fi — 4>i/y n'l and /2 = ^2/ y ?^2) 
and measure length in units of the coherence length ^1 = {2miVn[)~^^'^. Parameter r = 
W/ VUV is larger than 1, and k = miVn[/m2Un'2 is the ratio of two coherence lengths 
squared. One can solve these equations analytically in the limit k ^ 0, where two coherence 
lengths differ significantly or in the case r — > 00. In the limit k — > the first term in 



the second equation is negligiable, so that /2 = if rf{ > 1, and /2 = y 1 — r/f if rff < 
1. Substituting this solution into the first equation and using dfi/dx — F{f), one can 
reduce Eq.(16) to an integer able first order differential equation. The solution satisfying the 
boundary conditions, /i = 1 at x = 00 and /i = at x = — 00 is 

fi{x) = Q{x - Xi)tcinh{x/V2) + e{xi - x)——-—j- . -, (18) 

(r + 1)^'^ cosh[(r — l)^/^(a; — X2)\ 

( 2r 

Hx) = e{x, - 1^1 - + 1) ,,3h^[(, _ i)V2(, _ ,,)] j ' (19) 

where Q{x) is the ©-function, tanh(xi/-\/2) = r~^/^ and cosh[(r — iy^'^{xi — X2)] — 



2r/(r + 1). Interestingly the largest coherence length (in our case ^1) determines profile 
of both order parameters. The inter-component mutual repulsion (r) only slightly changes 
the characteristic lengths, as shown in Fig.2. 

The quasiparticle eigenstates of the phase segregated condensate are determined by the 
inhomogeneous BdG equations, Eqs.(3,4). Here we restrict the analysis by surface excita- 
tions at the boundary between two condensates in the limit r — > 00. Substituting Eqs. (18,19) 
into Eqs. (3,4) we obtain 4 coupled linear differential equations. In the limit r — > cxo the 
equations are decoupled in pairs, so that one can solve only first two equations at x > 0. 
Integration of the equations near zero provides the boundary condition rti(O) = fi(0) = 0. 
Then a simple analysis shows that there are no excitations localised near the boundary in 
that limit. 
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This conclusion seems to be quite general for any gapless Bose liquid. For example, 
one can consider a single component hard-core Bose gas with attractive impurity potential 
placed at the origin x = 0, so that the GP equation is 

+ / - + a6{x) = 0, (20) 

where o; > is the (dimensionless) strength of the potential. The solution is f{x) = 
coth[(|x| +xo)/V^] where sinh (2:0/2^/^) = 1/a. Substituting this solution into the BdG 
equations Eqs.(3,4) one can see that localised excitations do not appear even in that case. 
The condensate density increases at a; = and effectively screens out the attractive potential. 
However, localised excitations could appear in the gapped Bose-liquid, like the charged Bose 
gas. In that case they are formed below the plasma frequency. 

In conclusion, we have derived the Bogoliubov-de Gennes equations and the excitation 
spectrum of the two-component Bose-Einstein condensate with an arbitrary interaction be- 
tween bosons, including long-range and short range forces. Solving GP equations for segre- 
gated condensates we found the boundary profile, and showed that localised surface waves 
do not exist for strongly repulsive components. 

We greatly appreciate enlightening discussions of the nonlinear GP equations with 
Alexander Veselov. This work has been supported by EPSRC UK, grant R46977. 
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Figure captions. Fig. 1. The excitation energy spectrum of the two-component BEC 
with the long-range repulsive and hard-core interactions. In this plot rii — n2 — n, mi = m2 
and U = V = ujp, qi = q2 = e. Different curves correspond to the different values of 
W/V — 0.1 (dashed hne), 0.5 (dotted hne), and 0.95 (sohd line), respectively. Excitation 
energy is measured in units of the plasma frequency cUp — ATrne^/m and momentum k is 
measured in invers screening length Qs = (167re^nm)^/^. 

Fig. 2. Density profiles of two phase-separated condensates near the boundary for r = 2 
(solid hne) and r — 10 (dotted line), respectively. 



9 



